Spin Hall effect due to intersubband-induced spin-orbit interaction in symmetric quantum wells 
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We investigate the intrinsic spin Hall effect in two-dimensional electron gases in quantum wells with two 
subbands, where a new intersubband-induced spin-orbit coupling is operative. The bulk spin Hall conductivity 
cri, is calculated in the ballistic limit within the standard Kubo formalism in the presence of a magnetic field B 
and is found to remain finite in the B = limit, as long as only the lowest subband is occupied. Our calculated 
cjj, exhibits a non-monotonic behavior and can change its sign as the Fermi energy (the carrier areal density 
«2£>) is varied between the subband edges. We determine the magnitude of a* for realistic InSb quantum wells 
by performing a self-consistent calculation of the intersubband-induced spin-orbit coupling. 

PACS numbers: 72.25.Dc, 73.21.Fg, 71.70.Ej 



1. INTRODUCTION 

The Spin Hall effect (SHE)D refers to the spin accumula- 
tion with opposite polarizations at the two edges of a Hall 
bar, due to the transverse spin current induced by a driving 
longitudinal charge current. This effect was first theoretically 
proposed!] as arising from the spin-dependent scattering at im- 
purities with spin-orbit interaction. In the literature this is 
commonly referred to as the "extrinsic SHE" as it relies on the 
presence of impurities ("extrinsic mechanism"). In this case, 
the asymmetric Mott-skew and side jump scattering contribu- 
tions drive spin up and down electrons toward opposite direc- 
tions, thus giving rise to a net transverse spin current (with 
no charge current) and to spiruaccumulation at the edges. Re- 
cently, two theoretical workslill predicted that the spin-orbit 
effects on the band structure of semiconductors - the so called 
"intrinsic mechanism" for the SHE - can also give rise to a 
spin current perpendicular to an applied electric field, even 
in the absence of impurities. These authors have calculated 
the ballistic spin Hall conductivity (SHC) <x~ v , defined via 
J\ = ttj z J2 = h{v x cr z + o~ z v x )/4 = <j\ y Ey, where v x is the 
x component of the velocity operator and cr x>%z are the Pauli 
matrices, for a p-doped three-dimensional (bulk) valence band 
systemtl and for a two-dimensional electron gas (2DEG) with 
the Rashba spin orbit interactionU 

A number of theoretical papers have investigated the ro- 
bustness of the ballistic SHC cr| v as arising from the intrin- 
sic mechanism, against scattering by non-magneticlIDDDOD 
and magneticLJ impurities, its dependence on specific classes 
of SO interactions,!—! and th^ptetplay between the SO cou- 
pling .and. magnetic fieldsJiilOU Two early experimental 
effortsjjJiij have probed the spin Hall effect for electrons and 
holes.liZl Kato et cz/.L— I find that the SHE in n-GaAs epilayers 
is due to the extrinsic mechanism (non-magnetic impurities), 
while Wunderlich et aZ.l!j conclude that the SHE in the 2D 
hole gas probed in their experiment is intrinsic. Further in- 
vestigations have also shown that the SHE in 2DEGs is of the 
extrinsic type LI 

Following an early debate concerning the robustness of the 
intrinsic SHE, it is now well established that the dc SHC as de- 



fined above vanishes identically for model Hamiltonians with 
a linear-in-the-carrier-momentum SO interaction, such as that 
of Rashba and/or theJinearized Dresselhaus. This holds in 
both the ballistic caseLl ("clean limit") andmljbe.iiaut of weak 
scattering by non-magnetic impurities. UIHJjUlJLLI This re- 
sult can be understood by examining the relationship between 
the time derivative of the spin density and spin current in these 
systems. As pointed out in Refs. [| |], [Q, and [[]], for Rashba- 
type models da-^ldt oc f,, k — x,y and should vanish in the dc 
steady state regime (in the_nresence of some relaxation mech- 
anism) where dcr^/df = Oril 

The search for new materials that can exhibit the SHE as 
well as other types of SO interactions has continued over 
the yearshd Recently, we have introduced a new type of SO 
interaction present in III-V (or II- VI) zinc -blende -semicon- 
ductor quantum wells with more than one subbandLJU This 
intersubband-induced SO term is similar in form to the Rashba 
SO interaction. However, it couples electron states from dis- 
tinct subbands and hence can be non-zero even in structurally 
symmetric wells. For an electrpn-in a symmetric quantum 
well with two subbands we havetJtj 
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where m is the effective mass, <5 = {S e + fi )/2 and A£ = 
<5 - & e , with £> e and S a denoting the band edges of the lowest 
(even) and first excited (odd) subbands. t x >V:Z and cr xyz are 
Pauli matrices describing the subband degree of freedom and 
the electron spin, respectively. The intersubband-induced SO 
coupling 77 depends on the structural potential of the well, the 
electronic Hartree potential, and the external gate potential.U 
In this paper we calculate the dc spin Hall conductivity ojj v , 
for 2D electrons in the presence of the intersubband-induced 
SO interaction in wells with two subbands. We use the Kuba 
formula in the ballistic limit. We follow Rashba's approachU 
by performing our calculation in the presence of a perpen- 
dicular magnetic field B, which modifies the energy spectrum 
(Landau levels) thus allowing us to consistently include intra- 
and inter-branch transitions in the Kubo formula, and then tak- 
ing the B — > limit. With this procedure, known to produce 
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the correct vanishing of &*■ for the Rashba model,0 we de- 
rive an analytical expression for the dc SHC in two subband 
systems. More specifically, we find that: (i) cr z n , is non-zero 
for electrons whose Fermi energy Sf lies between the sub- 
band edges S e and S , i.e., when only the lowest subband is 
occupied, and (ii) cr z is null for Sf > S,„ i.e., when two sub- 
bands are occupied. Interestingly, the non-zero SHC is non- 
universal (e.g., it depends on 77), exhibits a non-monotonic be- 
havior and a sign change for S e < Sf < &o- In addition, the 
SHC presents finite jumps at the subband edges when plotted 
as a function of & F (or the carrier areal density «2£>X due to 
discontinuities in the density of states contributing to the spin 
density. We have also performed a detailed self-consistent cal- 
culation of the energy spectrum and wave functions for real- 
istic wells, from which we determine 77 and the corresponding 
a* . The magnitude of our self-consistently^determined SHC 
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is, however, much smaller than ctq = e/Az 

Before describing in detail our linear response calculation 
for the SHC, here we present a simple argument as to why 
we could expect a non-vanishing SHC for electrons in 2DEGs 
with two subbands where the intersubband SO interaction is 
operative. Despite the formal similarity between the SO term 
in our Hamiltonian, Eq. ([!]) and that of Rashba, here we find 
a different relationship between the spin current and spin den- 
sity, i.e., 
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in contrast to the Rashba model for which dcr x /dt oc f v In 
Eq. (^|) the spin current is related to the subband-related spin 
densities r A ® <t x and r v ® cr x . The second term in Eq. (Q), 
absent in the Rashba system, suggests that the (pseudo) spin 
density response to an applied electric field can contribute to 
the spin current - even in the dc steady state limit where d(r x ® 
o~ x )/dt = 0. Hence, 2DEGs formed in two subband wells 
with intersubband SO coupling may have a non-zero SHC. 
Our detailed linear response calculation below shows that this 
is indeed the case in the clean limit. 



2. MODEL HAMILTONIAN AND KUBO FORMULA 

Let us consider our Hamiltonian Eq. ([I]) in the presence 
of a magnetic field Bz = V x A by making the replace- 
ment p — > 7t = p - (e/c)A. For simplicity, we do not con- 
sider the Zeeman splitting (i.e, we assume a zero g factoio). 
Since our SO term couples only electrons with opposite 
spins in different subbands, the subband-spin Hilbert space 
{\bs); b — e,o,s =1, 1} can be divided into two independent 
subspaces T+ = {\e T>, \o |» and T- = {\o T>, \e |», and the 
2x2 Hamiltonian in each subspace A — + can be written as 



7t 



^ = — +6- 

2m 



AAS/2 T](n y +in x )/ti 
Tj(n y -in x )/n -AAS/2 



(3) 



The reduced Hamiltonian is identical to that of a 2DEG with 
a Rashba SO coupling of strength 7/ and an effective Zeeman 



splitting AAS. Note that only the sign of the effective Zeeman 
splitting differs between Ha m two subspaces. 

In the absence of the SO coupling (77 = 0), the Schrodinger 
equation gives rise to the Landau levels \bsn)o with the level 
index n > for each subband b and spin s. The Landau levels 
are evenly spaced by the cyclotron gap Tico c = \eB\jmc. The 
degeneracy (per unit volume) of each level and the magnetic 
length are, respectively, 1 /2nl 2 and / = ^chjeB. For non-zero 
SO coupling {r\ + 0) the coupling of Landau levels within the 
same subspace produces the mixed Landau levels \Ajj.n): 
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where (bi,b2) = (e, o) for A-+ and (o,e) for A--, with the 
corresponding eigenenergies 



&Ann = & + frUc (n ~ M 'An) ■ 



(5) 



Here ft — + denotes the spin branch. It is convenient to 
introduce dimensionless parameters S so = 2mrf/fi 3 u> c and 
S g = A&/hcj c for the SO coupling energy and the subband 
gap, respectively. In term of these we can define sin^,, = 
^nS S0 /y A „ and cos ^„ = £ a/7 An with & = (1 - AS g )/2 and 

J An = Jn&so + Note that the eigenstate for n = exists 
only for // = - sgn£). 

We determine the dc spin HalLconductivity at zero temper- 
ature by using the Kubo formulaU 
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where the primed sum indicates that it should be performed 
over the states with uiApn = £i/jn/^ < an d <^>An'n' > Sf/H. 
Here we have used the fact that the operators j\ and v v couple 
only states in same subspace and the matrices (Aiun\j\\Aiu'n') 
and (Aiu'n'\v y \Afm) are symmetric and antisymmetric, respec- 
tively. Since the matrix elements are independent of the guid- 
ing center position, the factor 1 /2nl 2 appears due to the Lan- 
dau level degeneracy. 

A more insightful analysis of the SHC can be achieved by 
expressing the operators in terms of commutators with the 
Hamiltonian. First, the relation v y = (i/h)[ r H,y] leads to 

(A/j.'n'\vy\Afxn) = H 2 {(i) All '„'-o)A l in) {An'n'\k x \Ann) , (7) 

with k x = (a 1 + a) I V2Z. This is a direct consequence of the 
definition of the ladder operator a — (I / ^f2Ti)(n x + m y ) = 

(y + Ppxlh + H 2 Py/K)/ V2/. From Eq. (§), on the other hand, 
we can obtain 
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+ 2^( Cl V' -w <V«') (^HTv ® o- x \Afi'n'} . (8) 
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Figure 1: Schematic diagrams describing the interbranch (long ar- 
rows) and intrabranch (short arrows) transitions. Left (a) and right 
(b) panels correspond to the cases £ F > fi„ and fi ( , < G F < B B , 
respectively. n+ denotes the highest filled Landau level in the spin 
branch pi. 



Substitution of Eqs. (|7|) and into Eq. (|6|) gives rise to (f- - 



r (t) 



+ 0-' 



(2) , with 



r (i) 



O"0 



r (2) 



=zz 

A imji'n' 



unrj 



(Apin\T x ® cr x \Api'n'} {Api'n'\k x \Apin) , 



(9a) 



tiAS (Afin\T y ® cr jr |/l y u'«') (/lyu'n'lfcjl/l/in) 



A fxnpi'n' 



mrj 



COAfin — l^A^'n' 



(9b) 

with <To s e/47r.D Each of the two terms can be further di- 
vided into two contributions coming from two different transi- 
tions classified by the restriction on the spin branches of states, 
to be summed over: interbranch and intrabranch transitions. U 
As can be seen from Fig. [|, the interbranch contributions 
describe transitions from the filled states in the lower branch 
pi — + to the empty ones in the upper branch pi = - and follow 
the selection rule \A+n) — » |A-n±l). The intrabranch tran- 
sitions are possible only in the vicinity of the Fermi level in 
each spin branch and obey a similar selection rule \Afin^ ll ) — > 
\Apirn ll +\), where n^, is the index of the highest filled Lan- 
dau level in the subspace A and spin branch pi; see Fig. [I]. In 
what follows, we write <x w = cr.£ + <x w , with i = 1 , 2 to 

' inter intra ' 

identify and separately investigate the contributions from the 
interbranch and intrabranch transitions. 



3. SPIN HALL CONDUCTIVITY 

First we investigate cr m . Making use of the fact that the 
matrices {Apin\T x ® cr x \Api'n') and {Api'n'\k x \Apin) are antisym- 
metric and symmetric, respectively, one can derive the two 
following identities 



where the asterisk over the sum indicates that it should be 
done over the states with S njl , 6^ < Sf- We have used 
the commutation relation [t x ® <x v , k x ] — to prove the sec- 
ond identity, which is valid for arbitrary A, pi, and n. Ap- 
plying these identities to Eq. (p>a|) reveals that <x (1) vanishes 
identically for arbitrary B and Sf ■ Interestingly, the inter- 
branch and intrabranch contributions exactly cancel out, i.e., 
(A , = -cr. i . This cancelation is similar to that which oc- 

mter intra r-i 

curs in the Rashba model.Lj Here, however, the extended op- 
erator t x ®ct x replaces the spin operator cr x in the Rashba case. 
More explicitly, we have 
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In the B — > limit, we obtain 



cr 



cro 



K\K2 + 1/2 



K\K2 
1 + 



K3 + Ki/2 



8 F >6 , 

Og<O^r<O 



(12) 



with K\ 
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A //c?/4 + K\Kr + 1/4. It is worth noting that cr ( : (or cr' 
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varies continuously as Sf passes through the upper subband 
energy edge S , even though the intrabranch transitions com- 
ing from the pi — - spin branch [see Fig. [jjb)] stop con- 
tributing at this energy. This is so because the contribution 
(nx-+\)ly, n + i in Eq. (11) vanishes as Sf — > 



We now evaluate the second term <x (2) of the spin Hall con- 
ductivity, Eq. (j9b|). In our model, <x (2) arises solely from the 
pseudo spin density response [see Eq. (|J)], which has no coun- 
terpart in the Rashba model. Differently than cr (1) , the expres- 
sion for cr (2) contains a factor a)^ n - <dah>„> in the denomi- 
nator [cf, Eqs. ( |9a[ ) an d (pbj )]; this prevents us from deriving 
identities such as Eq. (|l0[) for <x (2) . Hence, no exact cancela- 
tion between the interbranch and intrabranch contributions is 
guaranteed for <x (2) . In general, <x (2) is nonzero for arbitrary 
B. In addition, the denominator factor in <x (2) allows for the 
possibility of a finite contribution to <x (2) even if the matrix 
elements in the numerator of Eq. (|b|) vanish, provided that 
((^Afm ~ ^Afi'n') goes to zero as well. This, in contrast to the 
<x (1) case, makes the intrabranch contribution erf 2 ' discontin- 



{Apin\T x ®cr x \Api'n'){Api'n'\k x \Apin) = 0, (10a) uous as £ F crosses the subband edges. 



</IHTv ® cr x \Api'n') {Api'n'\k x \Apin) = 0, (10b) 



Explicitly, the interbranch and intrabranch contributions to 
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The intrabranch contrib ution from the fx = - spin branch, i.e., 
the second term in Eq. (13b), present only for Sf > S , does 
not vanish as Sf — > S and B — » 0. This leads to a discon- 
tinuity in cr (2) at £/r = S a . A similar disco ntinuity occurs at 
Sf 
&f 
to 



8 e due to the first term in Eq. (|13bJ) not vanishing as 
& e and B — > 0. In the B — > limit, Eq. (O) simplifies 
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Interestingly, the above results show that the interbranch 
and intrabranch contributions to cr (2) exactly cancel out in the 
B — > limit, provided that > £ , i.e., when both spin 
branches are filled. Note that this cancelation only occurs in 
the B — > limit. Hence, since cr (1) identically vanishes for any 
&f and B, we find that in the B — > limit crL is non-vanishing 
only for £ e < Sf <& being given by 



Hartree approximation - to determine the intersubband SO 
coupling strength rj for realistic symmetric wellsLI and the cor- 
responding spin Hall conductivity o%.(rj). We have considered 
a modulation-doped symmetric AlojIitojSb/InSb/AlojInojSb 
quantum well with only two subbands;LJ see Fig. g. The struc- 
ture comprises a single well of width L w = 15 nm and two 
n-doped semi-infinite adjacent regions. We assume that the 
ionized impurities giving up electrons to the well form a con- 
tinuum positive background with density p of width w = 5 nm 
(depletion layer). Our self consistent calculation follows that 
ofRef.||. 

In an earlier workJrJ we have found that narrow single 
InSb wells give sizable values for the intersubband-induced 
SO coupling strength rj. Here we self consistently calculate 
the corresponding-spin Hall conductivity as a function of the 
areal density «2D-LJ Figure || shows the calculated interbranch 



(center = (r. t + af{ ) and intrabranch {cr, 

\ miei lnter inter / v J 



intra intra 

conductivity contributions and the total spin Hall conductivity 
crL [see Eq. (JTsJ)] as a function of the areal density «2£> (or 
Sf) for a single InSb well. Note the two discontinuities of 
the SHC at densities corresponding to Sf = S (see vertical 
dashed line) and Sf = & e (at «2o = 0). The magnitude of 
these jumps are Acr/cr = atj /( 1 + ki) 2 and —K\j{\ - ki) 2 , re- 
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In the above expressions, the Fermi energy Sf (at B = 0) is 
given by 



(16) 



as can be straijghtforwardly derived from the B — spectrum 
of our system.U Next we calculate cr\ y for realistic wells. 
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4. SELF CONSISTENT CALCULATION AND RESULTS 

We have performed a detailed self-consistent calculation 
■ by solving the Schrodinger and Poisson's equations in the 



Figure 2: Self consistent potential profile of our InSb/AlojInn.ySb 
quantum well with two subbands. The calculated subband edges 
& e = 1.0 meV and £„ = 187.3 meV and the corresponding wave- 
funtions are also shown. The Fermi energy is G F = 100.3 meV. 
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Figure 3: Interbranch (dashed line) and intrabranch (long-dashed 
line) conductivity contributions and the total (solid line) spin Hall 
conductivity crL (in units of <tq = e/4n) as functions of the areal 
density w 2 d in the single InSb well shown in Fig. ^. Note the 
discontinuities at G F = fi„ and G F = G e , which correspond to 
ii2D = 10.3 x 10" cnT 2 and «2d = 0, respectively. 



spectively. These discontinuities come from the intrabranch 



contribution cr 

xy 

bution cr . , = 

xy .inter 



(A 2 }, , while the interbranch contri- 

lntra 

cr . varies smoothly with &f (recall that 



-(2) 



and 



= 0). In addition, the competition between or 
can lead to a sign change of crL as &f is varied, see 
Fig. Bl It should be noted that as the strength of the SO cou- 
pling j] (or K\) becomes smaller, Act and cri, diminish as well, 
thus vanishing completely at k\ = 0. 

Despite the interesting features displayed by the total spin 
Hall conductivity - the sign change and the non-monotonic 
behavior - as the electron density is varied between the sub- 
band edges, we find that the largest values of crL are ex- 
tremely small (~ 10~ 3 cro). 



[see Eq. (|9b|)]. The discontinuity Act in cr z xy can be explicitly 
related to the discontinuity Ap in the DOS as follows 
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Note that both expressions vanish as Ap/po — » 0, clearly 
showing that Act arises from the discontinuity in the DOS. 

Our intersubband-induced SO coupling mixes both the spin 
and the subband degrees of freedom simultaneously. There- 
fore it does not provide a mechanism to couple opposite spins 
within a given subband. This implies that any projection of 
our Hamiltonian into the lower subband, in the limit of large 
subband gap AS » <5 so , would not produce an effective SO 
coupling between spins within the lower subband so that no 
finite SHC appears. This limiting case study stresses that our 
finite SHC, even for Fermi energies near the lower subband 
bottom, is due to the subband-transfer process, even though it 
should be quite small in the limit AS » S so . We expect, how- 
ever, that additional subband-mixing but spin-preserving scat- 
tering mechanisms, such as impurity scattering, can induce 
an effective SO coupling within a given subband when medi- 
ated by the intersubband-induced SO coupling, which mixes 
both spins and subbands. This impurity-mediated SO cou- 
pling may then affect and possibly even enhance the strength 
of the SHC, as long as its momentum randomization is weak 
enough. 

As a final point we mention that interbranch contribution 
cr 1 inter reproduces the result for the SHC calculated via the 
Kubo formula in the absence of a magnetic field. Note that this 
quantity is non universal and vanishes as r\ — > 0. Nevertheless, 
we stress that it alone does not constitute the total SHC. 



6. SUMMARY AND FINAL REMARKS 



5. ADDITIONAL DISCUSSION 

Considering the formal similarity between our Hamilto- 
nian containing the intersubband-induced SO coupling and 
Rashba's, it is not entirely surprising that the SHC vanishes 
when both spin branches are filled. Novel effects of the 
intersubband-induced SO term arise when only the lower 
branch is filled: an incomplete cancelation between the contri- 
butions from the interbranch and intrabranch transitions lead 
to a nonzero SHC and discontinuities in it. The origin of 
these discontinuities can be traced back to the abrupt changes 
in B = density of states (DOS) of the model at the sub- 
band edges. The B — DOS for each spin branch is given 
by p T = |po(l + k\/2kt) with p = m/ntt 2 . Hence the 
DOS abruptly changes at &p = S a from a constant value 
p + + p_ = po to p+, thus giving rise to an abrupt loss of states 
contributing to the spin density response. In fact, the intra- 
branch contribution that is responsible for the discontinuity is 
proportional to the DOS, that is, 1/(0)^. - coahh^+i) k Pji 



We have calculated the ballistic spin Hall conductivity 
cr 1 for symmetric wells with two subbands in which the 
intersubband-induced SO interaction is present. We follow a 
linear response approach due to Rashba which consistently ac- 
counts for intrabranch and interbranch transitions in the Kubo 
formula. We find that crL is zero when the two subbands are 
occupied (similar to the Rashba model) and non-zero when 
only the lower subband is occupied. We have also performed 
a numerical self consistent calculation to determine the inter- 
subband SO strength for realistic InSb wells and have cal- 
culated the corresponding cr~ v . Even though the calculated 
cr^ shows interesting features such as discontinuities at the 
subband edges (due to discontinuities in the DOS), a non- 
monotonic behavior and a sign change as a function of the 
Fermi energy (or areal density), the magnitude of cr n , is much 
smaller than 1 in units of e/47r. 

It is conceivable that other materials systems, e.g. metal- 
lic surfaces and interfaces, can display intersubband induced 
SO interaction, in addition to the usual Rashba. For instance, 
metallic surfaces with unconventional spin topologyU where 
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deviations from the usual Rashba model have recently been 
reported (and for which the SO couplings are much stronger 
than in semiconductors) are a possibility. Perhaps in these 
systems the SHC due to the intersubband SO coupling would 
be sizable. 

One caveat of our calculation is that we use a spin current 
definition which is a simple extenskip for two subbands of the 
conventional (not-uniquely defined^]) symmetrized product of 
the spin and velocity operators as in Rashba model. Hence all 
the issues related to the reality of these currents and whether 
or not they would lead to spin accumulation in finite samples 
appear here as well More work is certainly needed to ad- 
dress these issues.L-1 Using non-equilibrium Green functions 
on a lattice, we have performed some simulations!-] of the 
spin density in bilayer systems with inter-layer SO orbit cou- 



pling, whose Hamiltonian maps onto our two-subband one. 
Our preliminary results show that the spin density changes as 
compared to the single layer case (Rashba model). Finally, we 
emphasize that the role of impurities, which we believe should 
not kill the effect discussed here, remains an interesting prob- 
lem for further investigations. 
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